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Abstract 

Measurements of heavy quark production and suppression in ultra-relativistic nuclear collisions probe the inter- 
actions of heavy quarks with the hot, dense medium created in the collisions. ATLAS has measured heavy quark 
production in yfsNN - 2.76 TeV Pb+Pb collisions via semi-leptonic decays of open heavy flavour hadrons to muons. 
Results are presented for the per-event muon yield as a function of muon transverse momentum, pj, over the range 
of 4 < < 14 GeV. The centrality dependence of the muon yields is characterized by the "central to peripheral" 
ratio, Rqp- Using this measure, muon production from heavy quark decays is found to be suppressed by a centrality- 
dependent factor that increases smoothly from peripheral to central collisions. Muon production is suppressed by 
approximately a factor of two in central collisions relative to peripheral collisions. Within the experimental errors, the 
observed suppression is independent of muon pj for all centraUties. 
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1. Introduction 

Collisions between lead ions at the LHC are thought to produce strongly interacting matter at temperatures well 
above the QCD critical temperature. At such temperatures, strongly interacting matter is expected to take the form of 
"quark-gluon plasma." High-px quarks and gluons generated in hard-scattering processes during the initial stages of 
the nuclear collisions are thought to lose energy in the quark-gluon plasma resulting in "jet quenching" 1 1 1 . Since 
the energy loss results from the interaction of a quark or gluon with the medium, jet quenching is thought to provide 
a valuable tool for probing the properties of the quark-gluon plasma L2j|3jl4jj. The contributions from radiative 15] 
and collisional ||6| energy loss in weakly coupled calculations are expected to depend on the heavy quark mass. In 
particular, the mass of heavy quarks is expected to reduce radiative energy loss for quark transverse momenta less than 
or comparable to the quark mass (;«), pj < m, through the dead-cone effect |7|. However, measurements of heavy 
quark production at RHIC via semi-leptonic decays to electrons showed a combined charm and bottom suppression 
in Au+Au collisions comparable to that observed for inclusive hadron production i[8ll9l [T0ll . There is disagreement in 
the theoretical literature regarding the interpretation of the RHIC heavy quark suppression measurements inTl[T2l[T3l 
particularly regarding the role of non-perturbative effects 1 14|. It is clear that measurements of heavy quark quenching 
at the LHC are an essential complement to inclusive jet or single hadron measurements. 

This paper presents results on heavy-quark production through muon semi-leptonic decay, using the ATLAS de- 
tector ifTSl and approximately 7 //b ' of Pb+Pb data at y/J^N - 2.76 TeV collected in 2010. Further details of the 
analysis can be found in 1 161 . The measurements were performed for several intervals of collision centrality over 
the muon transverse momentum range 4- < pj < 14- GeV and the yields compared to those in a peripheral bin using 
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Rrp, Rrp - .,'^°"ccni " — , where dn'^'^'" and dnP''"P'' represent the per-event differential rate for the same observable 

in central and peripheral collisions, respectively. (A^coii'^'^"') and (A^coii'"^"'''') represent (A^coii) values calculated for the 
corresponding centralities. 

2. Experimental setup and muon reconstruction 

Muons were detected by combining independent measurements of the muon trajectories from the inner detector 
(ID) and the muon spectrometer (MS), which cover \ri\ < 2.5 and \ri\ < 2.7 in pseudorapidity respectively. Two 
forward calorimeters placed symmetrically with respect to the collision point and covering 3.2 < \ri\ < 4.9 are used 
to characterize Pb+Pb collision centrality. Minimum bias Pb+Pb collisions were identified using measurements from 
the zero degree calorimeters (ZDCs) and the minimum-bias trigger scintillator (MBTS) counters. Cuts are imposed 
both online and offline to reject background and ensure selection of heavy ion collision events. In total, 53 million 
minimum-bias Pb-i-Pb events are selected. Previous studies |17| indicate that the criteria used to select minimum- 
bias hadronic Pb+Pb collisions have an efficiency of 98 ± 2%. The centrality of Pb+Pb collisions was characterized 
by I,E^'~^^, the total transverse energy measured in the forward calorimeters (FCal). A standard Glauber Monte- 
Carlo analysis was used to estimate the average number of participating nucleons, (A^pan), and the average number 
of nucleon-nucleon collisions, {Ncon), for Pb+Pb collisions in each of the centrality bins. The Rqp measurements 
presented in the paper use the 60-80% centrality bin as a common peripheral reference. 

The performance of the ATLAS detector and offline analysis in measuring muons in Pb+Pb collisions was eval- 
uated using a Monte Carlo (MC) data set |.18j . The MC data set was obtained by overlaying GEANT4-simulated 
lfT9ll - 2-76 TeV p+p dijet events on to I milUon GEANT4-simulated minimum-bias Pb+Pb events obtained from 
version 1.38b of the HIJING event generator EOl . HIJING was run with default parameters, except for the disabling 
of jet quenching. To simulate the effects of elliptic flow in Pb+Pb collisions, a parameterized cos 20 modulation 
that varies with centrality, and pj was imposed on the particles after generation. The efficiency for reconstructing 
muons associated with heavy flavour decays was evaluated using the MC sample described above for different bins in 
collision centrality. 



3. Muon signal extraction and results 

The measured muons consist of both "signal" muons and background muons. The signal muons are those that 
originate directly from the Pb+Pb collision, from vector meson decays, and from heavy quark decays. The back- 
ground muons arise from pion and kaon decays, muons produced in hadronic showers in the calorimeters, and mis- 
associations of MS and ID tracks. For this analysis a weighted combination of two discriminant quantities has been 
used to from a "composite" discriminant, C . These two quantities are: a) the fractional momentum imbalance, 
Apioss/piD, which quantifies the difference between the ID and MS measurements of the muon momentum after 
accounting for the energy loss of the muon in the calorimeters, and b) the "scattering significance", S , which charac- 
terizes deflections in the trajectory resulting from (e.g.) decays in flight. Further details are given in 1 16 1. 

The signal and background distribution of C were obtained separately from the MC and used to perform template 
fits to the C distribution of the data. Specifically, a probability density distribution for C, dP/dC, is formed assuming 

dP dP dP 
that a fraction, /s, of the measured muons is signal, — - /s— + (1 - /s) tt; , where df/dC|s and dP/dCle 

dC dC s dC B 

represent the Monte Carlo C distributions from signal and background sources, respectively. To account for possible 

differences in the C distribution between data and MC due to (e.g.) inaccuracies in the MC description of multiple 

scattering, the fits allow for limited adaptation of the MC templates to the data. In particular, the fit allows a shift 

and re-scaling of the C distribution, C = a + (C) + b{C - (C)), where (C) is the mean of the C distribution. In 

addition, a Gaussian smearing of the MC C distribution was included such that the data C distribution was fit to 



dF ^ / dP 
dC ~ dC 



-C'-/2iT- 

A kernel estimation method ll2T1l was used to produce the unbinned 



probability density distribution, dP' jdC from the MC signal and background samples. The signal fraction, /s, was 
then evaluated for each centrality and muon p-Y bin using unbinned maximum likelihood fits with four free parameters, 
/s, a, b, and cr. Examples of the resulting template fits are shown in the left panel of Fig.[T]for the 0-10% and 60-80% 
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Figure 1: Left: Examples of the template fits to measured C distributions (see text) in 0-10% (left) and 60-80% (right) centrality bins for two px 
intervals: 5 < px < 6 GeV (top) and 10 < px < 14 GeV (bottom). The curves show the contributions of signal and background sources based on 
the corresponding C distributions obtained from MC including the shift, re-scaling, and smearing modifications (see text). Right: Fractions of 
signal muons, fs , in the measured muon yields as a function of muon px in different bins of collision centrality. The points are shown at the mean 
transverse momentum of the muons in the given px bin. The shaded boxes indicate systematic errors, which can vary from point-to-point. The 
error bars show combined statistical and systematic uncertainties. Please refer to 1 16] for details. 

centrality bins and for two muon pj intervals. Typical fitted values for the a,b, and cr parameters are, a ~ 0.02, 
b ~ 0.95 - 1.05, and cr ~ 0.002. The combination of the MC signal and background C distributions are found to 
describe well the measured C distributions. The description remains good even if the adaptation of the template is not 
implemented, and the resultant variation in the fit contributes to the systematic error, as discussed below. 

The signal muon fractions extracted from the template fits for all pj and centrality bins are shown in the right 
panel of Fig. [T] The statistical uncertainties on /s from the fits represent Icr confidence intervals that account for the 
limited statistics in the data C distributions and correlations between the fit parameters. The uncertainties in the fit 
results due to the finite MC statistics were evaluated using a pseudo-experiment technique in which new pseudo-MC 
C distributions with the same number of counts as the original MC C distributions were obtained by statistically 
sampling the MC distributions. The resulting pseudo-MC distributions were then used to perform template fits. The 
procedure was repeated eight times for each pj and centrality bin and the standard deviation of the resulting /g values 
in each bin was combined in quadrature with the statistical uncertainty from the original fit to produce a combined 
statistical uncertainty on fs. The fractional uncertainties on fs due to MC statistics are found to be below 2% and are 
typically much smaller than the uncertainties from the template fits. 

The sensitivity of the measured fs values to the adaptation of the MC C distributions to the data was evaluated 
by performing the fits without implementing the adaption, i.e. by fixing the a and cr parameters to zero and b to one. 
To test the sensitivity to the primary hadron composition in the MC, and in particular to the relative proportion of 
kaons and pions (K/n), which may differ from that in the data, new MC background C distributions were obtained by 
separately doubling the n and K contributions. Potential systematic uncertainties resulting from the template fitting 
procedure were evaluated using a simple cut-based procedure applied to the C distributions. The total uncertainties 
from the above considerations are around 7% and 5% for 7 - 8 GeV and 10-14 GeV bins in the - 10 % centrality 
interval, and are around 19% and 6% for 7-8 GeV and 10 - 14 GeV bins in the 60 - 80 % centrality interval. 

The resulting central-to-peripheral values, Rqp, are shown in the left panel of Fig. |2] The Rcp values vary weakly 
with p-Y, and the points for each centrality interval are consistent with a -independent Rcp within the uncertainties 
on the points. The Rcp varies strongly with centrality, increasing from about 0.4 in the 0-10% centrality bin to about 
0.85 in the 40-60% bin. The variation of Rcp with centrality as characterized by the average number of participants. 
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Figure 2: Left: Muon Rcp as a function of pj for different centrality bins. The points are sliown at the mean transverse momentum of the muons 
in the given pj bin. The error bars include both statistical and systematic uncertainties. The contribution of the systematic uncertainties from S^oll 
and £, which are fully correlated between pj bins, are indicated by the shaded boxes. Right: Muon Rqp as a function of (A'part) for different bins 
in muon pj. The eiTor bars show combined statistical and systematic uncertainties. The sets of points for the different px bins are successively 
displaced by AA'pa,-t = 6 for clarity of presentation. Please refer to 1 1 6 1 for details. 



(A^part), is shown in the right panel of Fig. |2]for the different /?t bins included in the analysis. The Rqp decreases 
smoothly from peripheral to central collisions. The centrality dependence is observed to be the same for all pj bins 
within the experimental uncertainties. 

4. Conclusions 

This paper has presented ATLAS measurements of muon production and suppression in yfsNN - 2.76 TeV Pb-i-Pb 
collisions in a transverse momentum range dominated by heavy flavour decays, 4 < pj < 14 GeV, and over the 
pseudo-rapidity range \t]\ < 1.05. The fraction of prompt muons was estimated using template fits to the distribution 
of a quantity capable of distinguishing statistically between signal and background. The centrality dependence of 
muon production was characterized using the central-to-peripheral ratio, Rqp, calculated using the 60-80% centrality 
bin as a peripheral reference. The results for Rcp indicate a factor of about 2.5 suppression in the yield of muons in 
the most central (0-10%) collisions compared to the most peripheral collisions included in the analysis. No significant 
variation of Rcp with muon pj is observed. The Rcp decreases smoothly from peripheral to central collisions. 
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